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Abstract. Vortex dynamics in a multiband anisotropic superconducting such as the 
Fe based superconductors, is interesting and potentially important for applications. In 
this study we examine flux-creep data for fields along the second magnetization peak 
observed in M{H) curves of BaFe1.9Nio.1As2 for if ||c-axis, ||a6-planes and H forming 
a 45° angle with a5-planes. We find that the M-H loops taken from the different field 
directions can be collapsed onto a single universal curve at all temperatures with a 
simple scaling factor equivalent to the superconducting anisotropy, showing not only 
that the vortex pinning is isotropic, three dimensional and most likely related to point 
like defects. The resulting critical currents however, taken from the Bean model 
appears to show enhanced low field pinning for H\\c. The features in the vortex- 
dynamics also differ in different field orientations and show no direct correlation with 
the second magnetization peak Hp as is the case with a direct crossover in pinning 
regimes. Isofield plots of the scaled activation energy obtained from flux-creep data is 
found to be a smooth function of temperature as the Hp{T) line is crossed consistent 
with a single type of pinning regime operating at this field, independent of field 
orientation. The functional form of the Hp(T) lines in the resulting phase diagrams 
also support this view. 
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1. Introduction 

The study of vortex-dynamics in the novel pnictide superconductors [1] has attracted 
increasing attention [2l[3llll|5l[6l[71[8l|9l [101 E]) due to the considerable high- Tc of 
these compounds, when compared to the conventional superconductors, and because of 
similarities with the high- Tc cuprates superconductors [I2]. As in the cuprates[T3| HI], 
most of the pnictides systems exhibit the second magnetization peak, or fish-tail, in 
isothermic magnetization curves, as well as large fiux-creep, allowing the study of 
different regions of the vortex-phase diagram in certain detail. The second magnetization 
peak is associated with a maximum in the critical current when measured as a function 
of field at a fixed temperature. This phenomena is not yet completely understood in 
pnictides^. Although pnictides have similarities with the cuprates, such as the layered 
structure and antiferromagnetism of the precursor non-superconducting system [12], it 
is well established that superconductivity in pnictides has a multi-band character |1 5 [ [16] 
for which it is predicted the existence of non-Abrisosov vortices pT[ [T8l [T9] potentially 
leading to new effects in the vortex matter. In this work we study the vortex dynamics 
in a BaFe1.9Nio.1As2 superconductor single crystal, by means of isofield magnetization 
M[H) curves and magnetic relaxation M{t) curves obtained for iJ||c-axis, iJ||a6-planes 
and H forming a 45° angle with afe-planes. . 

2. Experimental 

Magnetization data were obtained by using commercial magnetometers: a 5T MPMS 
based on a superconducting quantum interference device (SQUID) was used for most of 
measurements with if ||a6-planes ; and a 9T PPMS was used for the other measurements 
including all data with ii||c-axis and and if-45°-a6-planes. The measurements were 
made after lowering the sample temperature from above T^. in zero applied magnetic field 
(ZFC-procedure). The studied sample, is a high-quality single crystal of BaFe1.9Nio.1As2 
with transition temperature T^. = 20 K, transition width ATc=0.3 K, mass = 120 mg 
and dimensions 0.5 x 1.6 x 0.02 cm. Details of the sample preparation can be found 
in Ref. [20]. For if||c-axis and H-45°-ab geometries, we carefully broke the sample 
and used a 43.1 mg piece with dimensions 0.55x0.5x0.02 cm. For if||a6-planes, the 
sample was attached to a hard-plastic slab perfectly inserted along the entire length of 
the straw tube used in the measurement systems, ensuring field alignment to within 2 
degrees with the a6-planes. 

3. Results and discussion 

Figure la shows a double plot of M{H) curves obtained at 15 K for H\\c and H\\ab 
where it is possible to observe the differences in the shape of the fish-tail between both 
field directions, which as shown below, are due to differences in the vortex dynamics. 
Figure lb shows M{H) curves at 16 K for all field directions, all exhibiting the second 
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Figure 1. a) double plot of M{H) at 15 K for H\\c and H\\ab.; b) M{H) curves at 16 
K for the three geometries. Insets: a) and upper b) shows a reduced plot of the same 
correspondent M{H) curves.; lower b), Jc{H) at 15 K for H\\c and H\\ab. Solid lines 
represent a fitting to the data. 

magnetization peak, where Hon represents the onset of the second magnetization peak 
with its maximum value occurring at Hp. The inset of Fig. la and the upper inset 
of Fig. lb show a reduced plot of the curves appearing in the respective main figures, 
where M is normalised by M{Hp) the value of magnetization at the respective Hp and 
equivalently the H scale by the respective irreversible field Hirr. These inset figures 
show that the M-H loops from the different field orientations can be collapsed onto one 
universal curve and therefore the same underlying physics is occurring independent of 
the orientation of field. If we compare Hirr{H\\ab) / Hirr{H\\c) , Hp{H\\ab) / Hp{H\\c) 
we find a factor that is of a value similar to the reported underlying anisotropy of the 
coherence length |27]. So the pinning properties are governed by the superconducting 
anisotropy. Although this appears straightforward, when we compare the absolute value 
of the critical current density from the different field directions (by using the Bean Model 
[2T| and the crystal dimensions above listed), as shown in the lower inset of Fig. lb for T 
= 15 K, we find that at low field, Jc{H\\c) exceeds Jc{H\\ab). Interestingly, this unusual 
effect in the critical current, is predicted in Ref. [H] for modest anisotropic systems 
with single pinning centers (acting as point-like pinning) which are physically large 
compared to the coherence length exhibiting strong pinning for fields applied parallel 
to the c-axis. This appears to be the case for low fields, implying that althought the 
field scales with the underlying anisotropy, the complexity of the multiband nature of 
the superconductor is evident in the absolute magnitude of the critical current density. 
In more anisotropic systems it is relevant to consider the 3D simulation of vortices 
interacting with single pinning centers for a fixed magnetic field presented in Refs. 
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[221 123], showing that an anisotropy dependent transition from 3D to a quasi 2D hne 
vortex appears associated with an increase in the critical current. Ahhough in the 
system we consider here with much more modest anisotropy a 3D to 2D crossover is 
not directly relevant. The lower inset of Fig. lb shows fittings of Jc{H) at 15 K to 
the expression J^{B)=A/[{B - Bpf + {l^Bff/^ of Ref. |30] where A, Bp and AS 
are fitting parameters, where Bp is the peak position and AB is the peak width. The 
resulting fittings conducted in a wide field range are excellent suggesting that the shape 
of the Jc{H) curve is primarily governed by vortex lattice softening as a function of field. 
The fitting shown in the lower inset of Fig. lb produced Bp = 11.7 kOe and AB = 20 
kOe for H\\c and Bp = 35.3 kOe and AB = 47 kOe. Similar fittings were obtained at 
different temperatures. 

The almost perfect symmetry of the M{H) curves of Fig. 1 with respect to the 
X-axis is evidence that bulk pinning is dominant. Also, the equilibrium magnetization, 
Meq, defined as the average value of M on both branches of a M{H) curve, is very 
small, so we may use values of M instead (M — Meq) in the analysis that follows. Near 
Tc for temperatures above 19 K, the second magnetization peak is no longer observed. 

In order to understand whether a vortex lattice softening model is sufficient to 
explain all features of the data, a more detailed vortex dynamics study was performed 
by collecting isofield magnetic relaxation data, M{t) curves, along several isothermic 
M[H) curves, for fields lying below and above the second magnetization peak, for the 
three geometries. We also obtained isofield magnetic relaxation curves as a function of 
temperature. Magnetic relaxation data were collect for 2 hrs when obtained in both 
branches of M{H) curves and for 3.5 hrs when only in the lower branch. We also 
measured long time magnetic relaxation for approximately 12 hrs in different regions of 
the M{H) curve for iJ||c-axis. All M{t) vs. log(t) curves show strictly linear behavior, 
starting above a transient time tq ~ 3-4 min. Such a large transient time was observed 
before in BaFei.82Nio.i8As2 with Tc=8 K, and seems to be intrinsic of the material [8]. 
The linear behavior with log(t) was also observed for the 12 hrs relaxation curves. This 
fact suggests that it is more appropriate to analyse the data by using the relaxation 
rate R = dM/dln{t) as defined in Ref. [21]. As shown in Ref. [21] one may obtain 
information on the apparent activation energy, but this quantity may not have any 
physical meaning in our data, since the magnetization Mq at the time t = 0, above 
which logarithmic relaxation should start, is not well defined, due to the ~ 4 min long 
transient region (see Fig. 2a below). 

Figure 2 shows results of the analysis of M{t) data for iJ||c-axis. The logarithmic 
relaxation of the magnetization, is exemplified by the selected curves plotted in Fig. 2a. 
Figure 2b shows a double plot of —R{H) vs. H and the corresponding M{H) curves 
where the 2hrs relaxation data are also plotted. In Fig. 2b, we plot -R as a function 
of increasing and decreasing field and we define Hi as the minimum in R [R increases 
above Hi) and H2 a maximum (i? decreases above H2). These two fields. Hi and H2, 
are defined from data obtained in the incresing field branch of M{H) curves. It is 
interesting to note that the resulting -R{H) curve resembles the correspondent M[H) 
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Figure 2. ||c-axis: a) selected magnetic relaxation data b)Double plots of ~R vs. H 
and M{t) vs. at T = 12 K 

curve after a shift to the right. The same trend as shown in Fig. 2b was observed for T 
= 14 K and 16 K data. 

Figure 3 shows the results of R vs. H for relaxation data obtained on the 
lower branch of all M{H) curves, where the arrows pointing up and down show the 
approximately positions of H2 and Hi respectively. The position of Hp in each curve of 
Fig. 3 (not shown) lies approximately at {H1+H2) /'i. It is clear from Figs. 2b and 3, 
which are representative data from all experiments, that there is no apparent change in 
the relaxation rate as the peak field Hp is crossed. For the sake of comparison, the upper 
inset of Fig. 3 shows double plots, as in Fig. 2b, of data obtained at the lower branch of 
an M{H) curve at 8K for H\\c (left inset) and data obtained in a YBaCuO sample (T^ 
~ 92 K) (right inset), with relaxation data obtained during 60 min (if||c-axis) at the 
lower branch of an M{H) curve at T = 60 K. A perfect matching between Hp and the 
field position of the maximum in R is clear for YBaCuO, which in this case represents 
a pinning crossover |25] taking place as Hp is crossed. A direct comparison of these two 
figures shown in the upper inset of Fig. 3, suggests that the second magnetization peak 
in these two systems arise from different mechanisms. 

Figure 4 shows the results for if||a6-planes. In that case R also decreases as field 
increases above Hon, reaching a minimum at Hi, which lies well below Hp. For each 
curve of Fig. 4, arrows pointing up shows the position of Hp and pointing down of Hi. 
A similar result was previously observed for an overdoped BaFei.82Nio.i8As2 with Tc = 
8 K for if||c-axis and if||a6-planes[8]. The upper inset of Fig. 4 shows a double plot of 
-R vs. H and the correspondent M{t, H) vs. H curve, exhbiting a kind of a shift to the 
right effect as observed in Fig. 2b for H\\c. It is interesting to note that the field Hi for 
H\\ab is more separated from the field Hon when compared to the if||c-axis case. To 
check whether this change in R only occur when if ||a6-planes, we measure the sample 
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Figure 3. iJ||c-axis: R vs. H for the lower branch of M{t,H) curves. Upper msets: 
left) Double plots of R vs. H and M{t) vs. H atT ^8 K. right) Double plot of R vs. H 
and M{t,H) for YBaCuO at T = 60 K. 



for an intermediate geometry, with H forming a 45° with a6-planes. The resuhs for T 
= 14 K is shown in the lower inset of Fig. 4 (the same trend is observed for T = 16 K), 
showing that the position of the field Hi relatively to Hon for this geometry is similar 
to that observed for if||a6-planes. One may associate the differences in the relaxation 
rate (as in the position of the field Hi relatively to Hon) as the field rotates away from 
the c-axis (where normalised relaxation rate 5* = {1/M)R shows that the creep rate is 
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Figure 4. R vs. H for the lower branch of M{t,H) curves for H\\ab-p\a.nes. Insets: 
upper) Double plots of —R vs. H and M(i, H) at T = 16 K with H\\ab\ lower) Double 
plots of R vs. H and M{t,H) at T = 14 K with HA5°ab. 



approximately three times faster for H\\ab, than H\\c), with the unexpected behavior of 
the critical current observed with the same rotation, as shown in the lower inset of Fig. 
lb, where Jc{H\\c) exceeds Jc{H\\ab at low fields). 

To check whether a pinning crossover, or a vortex phase transition, would perhaps 
become evident near Hp using a different approach, we obtained a systematic set of 
isofield magnetic relaxation data as a function of temperature for H\\c and H\\ab. Now, 
as in Ref. [26], we calculate the pinning activation energy for a set of isofield data M{t, T) 
using the expression, U{M)=-Tln{dM/dt)+CT where C=ln.{Bwa/7rd) is a constant, 
where B ^ H is the magnetic induction, u is an attempt frequency, a is the hop distance 
and d is the sample size. It is believed [26] that the isofield U{M,T)/ g{T/Tc) (where 
g{T/Tc) is an appropriated scaling function of ?7[26]) should be a smooth function of \M\ 
within a temperature region with same pinning mechanism. Figure 5 shows the results 
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Figure 5. U{M,T)/{1 - TjTcfl'^ vs. \M\ for H = 25 kOc, iJ||c-axis. The solid line 
represents a log\M\ curve. Inset: R vs. T for H = 20 and 25 kOe, i?||c-axis. 

of U{M,T)/g{T/Tc) plotted against \M\ for H\\c, as obtained for H = 25 kOe data 
with C = 10 and g{T/Tc) = (1 — T/Tc)^^'^. Similar values of the constant C were found 
for cuprates superconductors [26] and for pnictides|ll |8]. The almost perfect log|M| fit 
linking all the data in Fig. 5 suggests the existence of only one pinning mechanism 
over the entire AT range, for which Hp is located near 11 K. The same trend was 
observed for U{M,T)/{1 - T/T^f/^ vs. \M\ with H = 20 kOe and C = 10 for H\\c, 
and for H\\ab data with H = 15 kOe where instead g{T/Tc), we used q(T/T,v^) [26] with 
Tirr{15kOe)=19 K and C = 13 (not shown). The inset of Fig. 5a shows plots of R vs. T 
as obtained from isofield data with H = 20 and 25 kOe for H\\c, which do not show 
any visible effect near Hp located at 12 K and 11 K in each respective curve. One 
would expect some feature in the plot of R vs. T as Hp is crossed, either for a pinning 
crossover or for a vortex-lattice phase transition. The same behavior was observed on 
similar plots for if||a6-planes with if = 15 and 20 kOe. 

Figure 6 shows a Hvs.T diagram where characteristic fields of the three 
geometries are plotted against temperature after being divided (scaled) by 
^J {sin{6)y + {l/?>){cos{0)Y where is the angle between H and the a6-plane, and the 
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Figure 6. Hvs.T diagram for the three geometries after scale the Y-axis for H\\ab 
and H4:5°ab. Sohd hnes represents a fitting of the data. 

factor 3^Hp{H\\ab) / Hp{H\\c)^Hirr{H\\ab) / Hirr{H\\c) is of the order of the system 
(upper critical field) anisotropy [27]. The collapse of the Hon, Hirr and Hp lines is 
evident and suggests that the same underlying physics is occurring independent of the 
orientation of field. Each solid line in Fig. 6 follows a (1 — T/Tc)"^ dependence with 
m = 1.8 (Hon), 1.6 (Hirr), 1.4 (Hp), 1.7 {Hi for H\\ab and H45°ab), 1.2 {Hi H\\c) 
and 1.3 (-^2)- It is interesing to note the differences in the temperature behavior of 
the Hi lines in Fig. 6. Figure 6 also suggests a collapse of Hi values ( for H\\ab 
and H4:5°ab) which lie between Hon and Hp. The feature in R associated with the 
field Hi appears to be intrinsic of the studied system, since it was also observed in 
an over doped sample of the same system [8j. Values of H2 (for H\\c) lie between Hp 
and Hirr. One may associate the anisotropic vortex dynamics observed here with 
the anisotropic neutron spin resonance found on a similar optimally doped sample [28]. 
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Interestingly, isotropic vortex dynamics [H] and isotropic neutron spin resonance [2S] were 
found on overdoped samples of BaFe2-xNi2:As2. Although evidence suggests no crossover 
in pinning regime associated with Hp, the peak field position Hp is time dependent, 
because of the difference of relaxation rate as a function of field. The time dependence 
of Hp is also expected for crossover from elastic to plastic pinning regimes [25], but 
as expected the Hp line we measure here does not fit the predicted elastic to plastic 
crossover supporting our claim that the second magnetization peak in this system is 
unrelated to this particular mechanism. Nevertheless the fittings of Jc{H) shown in the 
inset of Fig. lb suggest that the second peak occurs due to softening of the vortex 
lattice occurring near Hp. This scenario, would accommodate the fact that the rate of 
relaxation increases with field above Hp (as shown in Figs. 2, 3 and 4 for all geometries) 
as is expected for example in the case of a disordered vortex phase [5]. In that case, R 
should eventually decay as the magnetic field approaches Hirr which explains why the 
field H2 occurs above Hp, differently than in YBaCuO (see the upper inset of Fig. 3). It 
should be mentioned that it seems that there are few reports in the literature on vortex 
phenomenology with respect to the system studied here, but experiments conduct on 
BaFe1.8Coo.2As2 indeed show a disordered vortex phase for intermediate-high fields [31] . 

4. Conclusions 

In conclusion, our study of vortex dynamics shows several interesting observations. 
The field scales governing the shape of the M{H) loops scale with the underlying 
superconducting anisotropy. The shape of the Jc{H) curve exhibiting the peak effect, can 
be described in terms of vortex lattice softening as a function of increased magnetic field. 
The differences in the relaxation rate with field orientation and the field anisotropy of 
Jc indicates that the pinning mechanism may be related to single ("point-like") pinning 
centers which are physically large compared to the coherence length [11] for H\\c, acting 
as strong pinning at low fields. No change in pinning regime is identified as a function 
of field near the second magnetization peak of M{H) curves in BaFe1.9Nio.1As2 either 
from the functional form of the Hp{T) curve or from the isofield activation energies. 

Acknowledgement s 

SSS, LG and ADA acknowledge support from CNPq and FAPERJ, LFC thanks the UK 
Funding Council the EPSRC grant EP/H040048. 

References 

[1] Kamihara Y, Watanabe T, Hirano M and Hosono H 2008 J. Am. Chem. Soc. 130 3296 
[2] Yang H, Ren C, Shan L and Wen H-H 2008 Phys. Rev. 78 092504 

[3] Moore J D, Cohen L F, Yeshurun Y, Caphn A D, Morrison K, Yates K A, McGilvery C M, Perkins 
J M, McComb D W, Trautmann C, Ren Z A, Yang J, Lu W., Dong X L and Zhao Z X 2009 
Supercond. Sci. Technol. 22 125023 



Vortex dynamics as a function of field orientation in BaFei,gNio,iAs2 



11 



[4] Salem-Sugui Jr. S, Ghivelder L, Alvarenga A D, Cohen L F, Yates K A, Morrison K, Pimentel Jr. 

J L, Luo H, Wang Z and Wen H H 2010 Phys. Rev. B 82 054513 
[5] Kopeliansky R, Shaulov A, Shapiro B Ya, Yeshurun Y, Rosenstein B, Tu J J, Li L J, Cao G H 

and Xu Z A 2010 Phys. Rev. B 81 092504 
[6] Pramanik A K, Harnagea L, Nacke C, Wolter A U B, Wurmehl S, Kataev V and Buchner B 2011 

Phys. Rev. 83 094502 
[7] Bhoi D, Mandal P, Choudhury P, Dash S and Banerjee A 2011 Phys. C 471 158 
[8] Salem-Sugui Jr. S, Ghivelder L, Alvarenga A D, Cohen L F, Luo H and Lu X 2011 Phys. Rev. B 

84 052510 

[9] Sharma S, Vinod K, Sundar C S and Bharathi A 2012 larXiv:1204.0339l 
[10] Bonura M, Giannini E, Viennois R and Senatore C 2012 Phys. Rev. B 85 134532 
[11] van der Beek C J, Konczykowski M and Prozorov R 2012 Supercond. Sci. Technol. 25 084010 
[12] Norman M R 2008 Physics 1 21 

[13] Yeshurun Y, Malozemoff A P, and Shaulov A 1996 Rev. Mod. Phys.68 911 

[14] Mcnon G I 2002 Phys. Rev. B 65 104527 

[15] Cvetkovic V and Tesanovic Z 2009 Europhys.Lett. 85 37002 

[16] Ding A, Richard P, Nakayama K, Sugawara K, Arakane T, Sekiba Y, Takayama A, Souma S, Sato 
T, Takahashi T, Wang Z, Dai X, Fang Z, Chen G F, Luo J L and Wang N L 2008 Europhys.Lett. 
8347001. 

[17] Babaev E 2002 Phys. Rev. Lett. 89 067001 

[18] Babaev E and Speight M 2005 Phys. Rev. B 72 180502 

[19] Cho Y M, and Zhang P M 2008 Eur. Phys. J. B 65 155 

[20] Chen Y, Lu X, Wang M, Luo H and Li S2011 Supercond. Sci. Technol. 24 065004 
[21] Bean C P 1964 Rev. Mod. Phys. 36 31 

[22] Olson C J, Zimanyi G T, Kolton A B and Gronbech- Jensen N 2000, Phys. Rev. Lett. 85 5416 
[23] Olson C J, Reichhardt C, Scalettar R T, Zimanyi G T and Gronbech- Jensen N 2003 Phys. C 384 
143 

[24] Beasley M R, Labush R and Weeb W W 1969 Phys. Rev. 181 682 

[25] Abulafia Y, Shaulov A, Wolfus Y, Prozorov R, Burlachkov L, Yeshurun Y, Majer D, Zeldov E, 
Wuhl H, Geshkenbein V B and Vinokur V M 1996 Phys. Rev. Lett. 77 1596 

[26] McHenry M E, Simizu S, Lessure H, Maley M P, Couher J Y, Tanaka I and Kojima H 1991 Phys. 
Rev. B 44 7614 

[27] Rey R I, Carballeira C, Mosqueira J, Salem-Sugui Jr. S, Alvarenga A D, Luo H, Lu X, Chen Y 

and Vidal F 2012 unpublished 
[28] Lipscombe O J, Harriger L W, Freeman P G, Enderle M, Zhang C, Wang M, Egami T, Hu J, 

Xiang T, Norman M R and Dai P 2010 Phys. Rev. B 82 064515 
[29] Liu M, Lester C, Kulda J, Lu X, Luo H, Wang M, Hayden S M and Dai P 2012 Phys. Rev. B 85 

214516 

[30] Rosenstein B, Shapiro B Ya, Shapiro I, Bruckental Y, Shaulov A and Yeshurun Y 2005 Phys. 
Rev.B 72 144512 

[31] Yin Y, Zech M, WiUiams T L , Wang X F, Wu G, Chen X H and Hoffman J E 2009 Phys. Rev. 
Lett. 102 097002 



